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Abstract. We consider scalar and spinorial perturbations on a background described by a 
2 = 3 three-dimensional Lifshitz black hole. We obtained the corresponding quasinormal modes 
which perfectly agree with the analytical result for the quasinormal frequency in the scalar case. 
The numerical results for the spinorial perturbations reinforce our conclusion on the stability 
of the model under these perturbations. We also calculate the area spectrum, which prove to 
be equally spaced, as an application of our results. 



1. Introduction 

New Massive Gravity (NMG) is a novel parity-preserving, unitary [Ij, power-counting super- 
renormalizable |2|, three-dimensional model describing the propagation of a massive positive- 
energy graviton with two polarization states of helicity ±2 in a Minkowski vacuum, whose 
linearized version is equivalent to the Pauli-Fierz theory for a massive spin-2 field in three 
dimensions. The action of NMG consists of a "wrong sign" Einstein-Hilbert term plus a quadratic 
curvature term given by a precise combination of the Ricci tensor and the curvature scalar, which 
introduces a mass parameter [3j. As with other massive gravity theories, NMG also admits black 
hole-type solutions with several asymptotics and additional parameters [HE]. Even though this 
last feature could challenge the usual Einstein-Hilbert gravity, it is seen that the definition of 
mass in this new type of black holes is a conserved charge computed from a combination of the 
black hole parameters, which satisfies the first law of thermodynamics. 

In this contribution we present the results of a study on the stability of the three-dimensional 
Lifshitz black hole found in the context of the so-called new massive gravity (NMG) [6] . These 
solutions exhibit the anisotropic scale invariance, t — )• \^t, x — )• Ax, where z is the dynamical 
critical exponent. Specifically, we deal with the solutions found for the particular case of z = 3 
and a precise value of the mass parameter. 

We aimed to determine the stability of this kind of black holes by calculating the quasinormal 
modes (QNM) and quasinormal frequencies (QNF) of scalar and spinorial matter fields in 
the probe limit, i.e, there are no backreaction effects upon the asymptotically Lifshitz black 
hole metric. QNM do not depend of the particular initial perturbation that excited the black 
holes, while QNF depend only on the parameters of the black hole. For this reason QNF are 
considered the "footprints" of this structure. The study of QNM has motivated the development 



of numerical and analytical techniques for their computation ([3 [U [9]). QNM are important 
in astrophysical grounds as a way of possible detection of black holes through the observation 
of its gravitational wave spectrum. They are also relevant in the context of gauge/gravity 
correspondence [10\ lllj . where the Hawking temperature of the black hole is related to the 
temperature of a thermal field theory defined at the boundary, and QNF can be related to 
the dual field theory relaxation time. In this work we focus in the application of QNM as 
probes of stability. Moreover, the QNF we obtain are useful to calculate the area spectrum of 
Lifshitz black holes. Actually, according to Bekenstein the horizon area of a black hole must 
be quantized jI2J . so that the area spectrum has the form An = jnh, with 7 a dimensionless 
constant to be determined. In this way, we find this spectrum by using an approach suggested 
by Maggiore |I3j . who claims that QNM should be described as damped harmonic oscillators, 
and the module of the entire QNF is the proper physical frequency that should be used in the 
Bohr-Sommerfeld quantization of the adiabatic invariant I = f dE/uj[E). 

In the following sections we describe the model we dealt with, the equations for the scalar 
and spinorial perturbations, the numerical analysis, and the calculation of the area spectrum. 
In the last section we discuss our results and conclude. 



2. The black hole solution 

The black hole we study along this paper is a solution of the NMG (2 + l)-dimensional action. 



S 



1 



-9 



R-2X 



1 



ii-^un, — -K 



(1) 



where m is the so-called "relative" mass parameter, and A is the three-dimensional cosmological 
constant. This theory has a black hole solution with Lifshitz scaling z = 3 given by 
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with M an integration constant and = = ~2m^" ^1^°' NMG admits as a solution, the 
well-known Bahados-Teitelboim-Zanelli (BTZ) black hole with the dynamical critical exponent 
z = 1. The Kretschmann scalar for the metric ([2]), 
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diverges for r — ?■ 0. Thus, the black hole solution has a genuine spacetime singularity at the 
origin. From the Penrose-Carter diagram in Fig(T]we see that the spacetime singularity at r = 
is light-like covered by a regular event horizon at r = r-(_ = 



3. Scalar perturbation equation 

In order to consider the response of the Lifshitz black hole background to a scalar perturbation 
we consider a scalar field obeying the Klein-Gordon equation, 

□ $ = -LSm (V^/'^5,v) $ = , (6) 



Figure 1. Penrose-Carter diagram for the Lifshitz black hole. The singularity at r = is 
light- like and covered by a regular event horizon r_|_. 



where m is the mass of the field r, </>). After decomposing the angular part of the field this 
equation takes the following form, 

+ = 0. (7) 

Although this equation has an analytical solution, which we will show later, it is a useful test 
for our numerical code to check the numerical results. With this goal we need to rewrite ([7]) in 
a Shrodinger-like form, 

-d^X + dlX = V{r)X, (8) 



where is the so-called tortoise coordinate given by 

r* = I 
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and the scalar effective potential can be written as 

Now let us return to Eq.© and discuss its analytical solution. First of all we set the time 
dependence of the scalar field as ^'(t,r) = i?(r)e~*'^*, and then redefine the radial coordinate as 
r = l\fM jy. The resultant equation has the following solution in terms of the Heun confluent 



functions, 

+C,y'--{l-yr/'Rennc(^0,-a,(5,-^,^ + ^y^ , (11) 

where Ci and C2 are integration constants, while a = -v/4 + m?P^ and j3 = —iluj/M'^/'^. After 
imposing Dirichlet condition at infinity and the boundary condition of ingoing waves at the 
horizon we obtain the quasinormal frequency spectrum given by 
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where is a positive integer. These same frequencies were also obtained in subsequent 
works [nillSj by using equivalent analytical procedures. The imaginary part of the fundamental 
frequency [N = 0) is negative provided that 



y 7 + -m2/2 + — + 3v/4 + m2/2 > 1 + ^4 + ^2;2 _ (13) 
While the asymptotic frequency (A^ — )■ 00) is given by 

M3/2 

cj^ = -2(^/6 - 2) i — ^Af < . (14) 

Thus, since the imaginary part of the quasinormal frequencies is negative provided that the 
parameters respect the relation (jl3p . we can conclude that the model is stable under scalar 
perturbations. 

4. Spinorial perturbation equation 

Let us consider a spinorial field perturbing the Lifshitz black hole background, which obeys the 
covariant Dirac equation, 

i7We(^)''V^^-;U,^ = 0, (15) 

where Greek indices refer to spacetime coordinates {t,r,(p), and the Latin indices enclosed in 
parentheses describe the flat tangent space in which the triad basis e^^^^ is defined. The spinor 
covariant derivative is given by 

= 9^ + ^a;^('^)W [7(a), »] , (16) 
where uj^^""^^^^ is the spin connection, 

(a)(6) ^ ^^ia)Q^^ib)u ^ e^(a)r-^^e-(^') , (17) 

T'^fj_„ are the metric connections, and the gamma matrices we employ along this work are 
= i(72, 7^^^ = and 7^2) = g-^_ 



In order to solve the Dirac equation in a numerical way once again we should adapt it to a 
Shrodinger-like form. First, we redefine the two component spinor as, 



as 
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And in addition we make X± = Rj^ it to finally have 

(dl+u:^)x± = v±x±, 

where V± are the superpartner potentials, 

n dW 

written in terms of the so-called superpotential VK, 
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In the case of a massless spinor {pis = 0) the superpartner potentials reduce to 
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5. Numerical results 

Although in the scalar case we already found the corresponding QNF, the numerical analysis 
is an interesting tool to verify the applicability of certain numerical methods in asymptotically 
Lifshitz spacetimes. Specifically, here we use the Finite Difference and the Horowitz-Hubeny |16] 
methods. Let us begin with the scalar perturbation. By rewriting Eq. pup in terms of a new 
variable z = r"^, we obtain 
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where Zh = rf^- The parable in brackets tends to infinity as long as ^| + rn^P^ > 0, which is 
consistent with the Breitenlohner-Freedman-type bound for the present case. The roots of this 
polynomial potential are given by 



zo = 0, 

z+ = Zh, (26) 




Analytic 
Numeric 
Horowltz-Hubeny 



1.1 1.2 1.3 1.4 1.5 1.6 



Figure 2. Decay of scalar field with mass Figure 3. Imaginary part of scalar 
m = 1 and / = 1 for different values of black quasinormal frequencies. We display the 
hole mass M. results obtained using different approaches. 



If m?l'^ > —1, we see that Z- < 2;+. Thus, going back to the original variable r, the roots of 
the potential are r = with double multiplicity, r = -^/zT and r = (r = —^/zZ and r = —r^ 
are excluded as r > 0). Then, since is the biggest root and V{r') tends to 00 when r tends 
to 00, the potential is positive-definite in the region (r/i,oo). Therefore, the quasinormal modes 
for rrP'T?' > — 1 are necessarily stable [16]. 

In Figsi2]l3]we show the numerical results for the decay of the scalar field and a comparison of 
the scalar QNF obtained by three different approaches, analytical calculation, finite difference, 
and Horowitz-Hubeny methods. Our results reinforce the conclusion already found analytically, 
the z = 3 Lifshitz black hole is stable under scalar perturbations. Moreover, according to Figj3l 
the numerical results have a very good agreement with the analytical calculation. However, 
the Horowitz-Hubeny method gives unreliable results. We discarded explanations related to the 
not-convergence of the frequencies or the size of the black holes ocurring in other cases ]yi\ [18] , 
in spite of that, there is no clear explanation for the limitation of the method. In our case 
the asymptotic behavior of the black hole under study might play an important role in the 
convergence of the method. Nonetheless, a general criteria for the convergence of the Horowitz- 
Hubeny method remains an open question. 

Now turning to the case of spinorial perturbations, we can analogously analyze the 
superpartner potentials (|24|) and its derivatives, which we rewrite as 
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00 if ml > 0. 



We see that outside the event horizon V+ is positive-definite if ml > 0, and lim. 
if ml < 0. Whereas V- is positive-definite if ml < 0, and linir^ooV-i''') 
Furthermore, we notice that if ml = 0, we have a free-particle case. The decaying behavior of 
the massless spinor is given in FigsHHS] Thus, we conclude that the z = 3 Lifshitz black hole is 
stable under massless spinorial perturbations. 



Figure 4. Decay of massless spinor with 
/ = 1 and black hole mass M = 1.0 for 
different values of the azimuthal parameter 
m. 
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Figure 5. Decay of massless spinor with 
/ = 1 and black hole mass M = 1.5 for 
different values of the azimuthal parameter 
m. 



6. Applying QNM to thermodynamics 

In this section we use our analytic results for the QNF to find the area spectrum of the Lifshitz 
black hole. In order to apply Maggiore's method we need the asymptotic value of the QNF given 
in (fH|) . so that the proper physical frequency of the damped harmonic oscillator equivalent to 
the black hole QNM turns out to be 

Wp = 2(V6 - 2)— ^iV. (30) 

Now we can calculate the adiabatic invariant 
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where M = M'^ /2 is the ADM mass of the Lifshitz black hole, and Aw is the change of proper 
frequency between two neighboring modes, which can be obtained from (|30p . Thus, 
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which is quantized via Bohr-Sommerfeld quantization in the semiclassical limit. Then, expressing 
this equation in terms of the horizon area A = 27rr+, with r+ = we finally obtain the 

spectrum we look for, 

A = 27r{V6 - 2)nn , (33) 

with n an integer number. As we can see, the horizon area of the Lifshitz black hole is 
quantized and equally spaced, a result not expected from a theory containing higher order 
curvature corrections. In addition, since this black hole has an entropy proportional to its 
horizon area \19\ I20j . our result (j33p also says that the entropy is quantized with an evenly 
spaced spectrum. Subsequent studies on area spectrum of other kind of Lifshitz black holes [H] 
seem to confirm this last conclusion. 



7. Conclusions 

We have studied the stability of the three-dimensional Lifshitz black hole under scalar and 
spinorial perturbations in the probe limit through the computation of quasinormal modes. In 
addition, we have found the event horizon area quantization as an application of the results for 
quasinormal modes using Maggiore's prescription. 

Regarding the scalar perturbations, we found that the QNF are purely imaginary and negative 
pointing out a stable equilibrium configuration for the Lifshitz black hole. Our numerical results 
perfectly agree with the analytical calculation in the case of the finite difference method, which 
allows to obtain the temporal behavior of the fields showing all the stages of the decay. However, 
in the case of the Horowitz-Hub eny method the frequencies values fail to be obtained. An 
important observation is that apart from the numerical factor, the asymptotic scalar frequency 
is the same as that calculated in the hydrodynamic limit in the context of Gauge/Gravity 
duality [2T] . 

Concerning the spinorial perturbation, our numerical results also show probe massless spinors 
decaying, and therefore, the z = 3 Lifshitz black hole is stable under spinorial perturbations too. 

Additionally, using the asymptotical form of the analytical QNF we calculated the area 
spectrum of the Lifshitz black hole using Maggiore's method. Our results show that the horizon 
area and the entropy are quantized and their spectra are equally spaced. 

At last, we should stress that the definite answer on the stability of the Lifshitz black 
hole should come from gravitational perturbations, mainly from the tensor part of the metric 
perturbation, since NMG admits the propagation of gravitational waves contrary to Einstein 
gravity in three dimensions, that has no propagating degrees of freedom. 
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